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Abstract. We consider the one dimensional classical Ising model in a symmetric 
dichotomous random field. The problem is reduced to a random iterated function 
system for an effective field. The Dq-spectrum of the invariant measure of this 
effective field exhibits a sharp drop of all D q with q < at some critical strength 
of the random field. We introduce the concept of orbits which naturally group the 
points of the support of the invariant measure. We then show that the pointwise 
dimension at all points of an orbit has the same value and calculate it for a 
class of periodic orbits and their so-called offshoots as well as for generic orbits 
in the non-overlapping case. The sharp drop in the D 9 -spectrum is analytically 
explained by a drastic change of the scaling properties of the measure near the 
points of a certain periodic orbit at a critical strength of the random field which 
is explicitly given. A similar drastic change near the points of a special family 
of periodic orbits explains a second, hitherto unnoticed transition in the D q - 
spectrum. As it turns out, a decisive role in this mechanism is played by a specific 
offshoot. We furthermore give rigorous upper and/or lower bounds on all D q in a 
wide parameter range. In most cases the numerically obtained D q coincide with 
either the upper or the lower bound. The results in this paper are relevant for 
the understanding of random iterated function systems in the case of moderate 
overlap in which periodic orbits with weak singularity can play a decisive role. 



PACS numbers: 05.45.Df, 05.50.+q, 75.10.Nr, 05.70.Fh 



1. Introduction 

The properties of multifractal measures have attracted a lot of interest over the past 
two decades. Multifractals naturally appear in a variety of physical and mathematical 
contexts. From the beginning the one dimensional random field JjJ-flTj and random 
exchange [p^|-[pp| Ising models were prominent examples. In treating these systems 
some reduction scheme for the partition function like the transfer matrix method 
|Q, ||, [I?], |l8| or a method introduced by Rujan ^ usually is used [f|-(l^] to obtain 
a random iterated function system (RIFS) for a local effective field. This leads via 
a Frobenius-Perron or Chapman-Kolmogorov equation to an invariant measure which 
typically is a multifractal. Similar structures arise in other one dimensional disordered 
systems like phonons (2l[] or electrons |^2|, in random potentials, cf. also ]2^ j. 

In the early investigations of the random field Ising model the language of 
multifractals had not yet been developed and the results focused apart from calculating 
the free energy on the structure of the support of the invariant measure and the 

ground state properties of the system j|, ^, [| [l^, |l8| . More recently the uniqueness of 



Orbits and phase transitions in the multifractal spectrum 



2 



Gibbs measures and exact ground state properties were investigated p6[ . A connection 
to domain theory was established in |2f| . 

The invariant measure of the local effective field in the random field Ising model 
is a multifractal In various works the generalized box dimensions (generalized 
Renyi dimensions) D q p6{ of this measure were calculated for special q 4 , M , with 



perturbation expansions"^, [Tl], ^7) or numerical approximations jlO[ [l3, 27fl. Other 
authors focused on different concepts like the order-g free energy and its fluctuations 
p9[ p0[ or correlation functions Q . 

The systematic numerical investigation of the dependence of D q on the strength 
of the local random fields jl3) revealed the surprising feature of discontinuities (phase 
transitions) in the D q with negative q. 

Almost all features of the _D 9 -spectrum have by now been understood analytically. 
The most drastic transition in the -Dq-spectrum, the sharp drop of all generalized 

(2) 

dimensions D q with q < at some critical field strength he , also present in the context 
of a special model of neural networks |28L |2j| , was explained on a phenomenological 
level by the disappearance of deep cuts in the measure density at h < h ( c 2) @, |§. 
This disappearance of deep cuts in the measure density can be explained analytically 
by close investigation of the obtained nonlinear RIFS. 



In this paper we complete the analysis of the transition begun in llq, 30 and 



(2) 

explicitly calculate the critical field strength h c of the transition. The result is 
obtained by generalizing the analysis of the singularity (pointwise dimension) at fixed 
points to the singularity of orbits. Further application to a special family of periodic 
orbits explains a so far unnoticed smaller drop in the D g -spectrum at a critical 

field strength h^ a ^ which became observable because of increased precision in the 
numerical generation of the -D^-spectrum. Furthermore, the concept of orbits and 
their singularity also allows to give bounds on D q for any q and the exact value of 
D± oa in a wide parameter region of the random field strength h. The computation of 
D±oo generalizes earlier results in [fL0|| . 

Similar approaches and arguments may be found in the mathematical literature. 
In |3l) parabolic function systems with overlaps are considered, |}2j concentrates on 
measures obtained by infinite Bernoulli convolutions and |]33f investigates generalized 
dimensions D q of measures on general self-affine sets. 

In the following we consider the one-dimensional random field Ising model p|- p7t 
with the Hamiltonian 

AT-l N 
H N = - J ) SjSi+l - ^ h i S i' W 
i=l i=l 

in which Si denotes the classical spin at site i which takes values I or —1 and J is 
the exchange energy of adjacent spins. The local magnetic fields {hi} at the sites 
i = 1,...,N are independent identically distributed random variables. We restrict 
ourselves to dichotomous symmetric distributions, i.e. to probability densities with 
Dirac masses at ±/i, 

p(h i ) = ^5(h i -h) + ^5(h i + h), heR+. (2) 

An iterative reformulation of the canonical partition function yields the partition 
function of a single spin in an effective external random field x[ N ^ (the effective field 
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at site 1 in a chain of N spins) which is given by an iterative map Q , 

xW = hi+A(x%l), ^ = (3) 
1 /cosh/3(n: + J) 
(X) ~ 2j3 \ cosh 0(x — J) 

where (3 denotes the inverse temperaturej| The iteration is illustrated in figure |l|. As 
a shorthand we introduce 



f a (x) :=(Th + A(x), ere {+,-}, (4) 

such that the recursion (^) reads x^ — /o-j^+i) with hi —: Uih. By the 
reformulation of the canonical partition function we are thus led to a RIFS with 
smooth, strictly monotonously growing, contractive functions {/+, /-}= {A+h, A— h} 
and probabilities {p+,p-} = ^}. When viewing (Q) as a RIFS we will also write 

x n instead of x\ N ^ for the value of the effective field after n = N — i + 1 iterations. 
Please note that the transition from N to N + 1 spins implies prepending functions 
to the composition of functions, i.e. we need to consider x„ = f ai o . . . o fa- n (x), 

Xn+l = /<ri ° ••• ° fcr n ° f<r n+1 (x), etc. 

We introduce a symbolic dynamic in the obvious way: Let £„ be the set of finite 
sequences {cr}„ of n symbols Oi G {+, — }, i = 1, . . . ,n and Eoo the set of all infinite 
sequences {<r}. Given {a} we will write {o~} n for the head of the n leftmost symbols 
in {<r}. By we denote the composition of the n functions f ai , i — 1, . . . , n, i.e. 

ficr} n = firi ° fa 2 ° • • • /<t„ • The above mentioned properties of /+ and /_ imply the 
following facts: 

• The fixed points x* + and x*_ with f + {x* + ) = x* + and f^(x*_) = x*_ exist. 

• The interval / = [x*_, x*^} is the smallest interval with f a {I) C I. 

• The limit linin^oc f{ a } n {xQ) exists for any {a} € E^ and xo <E I and does not 
depend on xq. We thus can define x*^ := lim^oc f{cr} n {xo) and the (constant) 
function /{ CT }(a;) := a^i.- Through this definition air , is the unique 
fixed point of f{ a }- 

• We denote the fixed points of finite compositions by x*^ . If {a} is periodic 
with period n then Xr CT | = x \ a ] ■ 

• We name the n-fold images I{ a } n '■— f{a} n {I) °f the invariant interval / bands 
(of order n). For given {a} they have the inclusion property C for 
n > m. Furthermore, the construction implies £/ CT | £ -^{>}„ f°r au n - 

• If the first order bands i+ and J_ do not overlap, 1+ n /_ = 0, none of the higher 
order bands overlap, Ij a } n HZrga = We call this the non-overlapping case and 
denote the first order gap by A := [/_(x+), f+(x*_)]. The inverse statement is 
also true. If the first order bands overlap the corresponding higher order bands 
also overlap. This situation is illustrated in figure [jj. We denote the first order 
overlap by O := I + n J_ . 

The RIFS (plpl) induces a probability density p n for the effective field x n . We 
write P n {x) := Jq Pn(€)d£ for the corresponding distribution function and fj. n for 

i Note, that we find it convenient to use a slightly different notation than in previous work as e.g. 
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Figure 1. (a) Mapping in the case of overlapping bands 7+ = /+(/) and 
/_ = /_ (/). In (b) the first few images of the interval I are shown. The increasing 
complexity of the band structure is obvious (h = 0.55, (3 = 1, J = 1). 



the corresponding measure. P n , p n and \i n can iteratively be constructed using the 
Frobenius-Perron (Chapman-Kolmogorov) equation induced by (|^) which reads for P„ 



Pn(x) 



f dh p(h)P n ^(A-\ X - *0) = £ \Pn-l(f-\x)) 



(•5) 



with Po(x) = Q(x), being the Heaviside function. This choice of Pq encodes the 
free boundary conditions chosen in (^) . Again the properties of /+ and /_ imply some 
direct consequences: 

• The Frobenius-Perron equation has a unique fixed point Poo (/too) Hj @ ano - the 
invariant measure /i^ is ergodic. 

• The reiterated application of the Frobenius-Perron equation to an arbitrary initial 
measure /iq (distribution Poo) converges to //oo (Pn) in Hutchinson topology or, 
as / is compact, equivalently in the weak topology of measures p4j. Thus, /x^o 
(Poo) is the measure (distribution) of the effective field x in the thermodynamic 
limit N — > oo for any boundary condition. 

• The explicit form of the n-th iterate P n (x) reads 



Pn(x) 



£ 



(6) 



which is in the limit n — ► oo a path integral in the space of symbolic dynamics. In 
the non-overlapping case the sum on the right hand side of (^|) has only one term 
for each x. In the overlapping case typically more than one term contributes for 
each x . There are however x € supp /i^o for which still only one term contributes. 

The support supp fi^ C I is the attractor of the PJFS {/+,/-} J34|]. 

For any x € supp /loo one can find a sequence {a} with x — %{ a y In the non- 
overlapping case this relation between and supp fioo is one to one. 
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In the following we closely investigate the multifractal properties of the invariant 
measure /ioo- To this end we study the generalized box counting dimensions 

D q = J- lim l^il, (7) 
q q - 1 e^o In e 

where fit = fJ.oa(B E (xi)) are the measures of boxes (intervals) of size e covering 
supp /ioo , and the pointwise dimension 



at individual points x € supp /ioo. We will synonymously use the singularity a := 
D p — 1. The intricate interplay between these quantities will explain the transitions in 
the -D ? -spectrum mentioned above. The obtained characterization of /ioo the measure 
of the effective field x, is a necessary prerequisite for the more complicated treatment 
of the distribution of physical quantities like the local magnetisation which will follow 
in a later publication. 



The paper is organised as follows. In subsection LJ we briefly review the known 
results about prominent features of the -Dg-spectrum. In section || the concept of 



orbits and their singularity is introduced and in 2.1 this singularity is calculated for a 



class of periodic orbits and their so-called offshoots. In subsection |2.2| we calculate the 



generic singularity of arbitrary orbits in the non-overlapping case. In subsection 2.3 



2) 



we treat the overlapping case and discuss effects occurring if points of an orbit enter 
the overlap. The results are used to explain the transition in the -Dg-spectrum at h 

(2) 

and calculate he 'as a function of temperature T and coupling strength J explicitly. 



We then apply a similar analysis to the transition at h^c°^ in subsection |2.4| . In section 
[| we give the extended lower and upper bounds on the -Dg-spectrum obtained from 
the analysis of the singularity of specific orbits. Finally, some conclusions are drawn 
in section [|. 

1.1. Known results 

In this subsection we summarize in short previous work on phase transitions in the 
invariant measure /ioo- For large h the support of /ioo is non-connected and similar to 
a multi-scale Cantor set 



At a critical value he of h the support of /ioo becomes connected for all h < h c 
The value of h^ is determined by the overlap condition for the first bands, 
f-(x* + ) = f+{x*_). This results in [|| 

^^arcosh^-l)/^. (9) 

The transition can be seen in the densities p n of the approximations fi n of the invariant 
measure /ioo | fl3l . In the -Dg-spectrum the transition is visible as the point where Dq 
becomes 1 , cf . figure ||. 

At h^ < h{ lS> the invariant measure density jumps from infinity to zero at the 
boundary x* + and x*_ of / |l3| . This is an effect solely depending on the scaling of 
the measure at the fixed points x*_ and x* + . The measure density p n (x±) diverges for 
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Figure 2. 

local effective field with q = -20, -6, -3, -2, -1, 0, 1, 2, 4, 20 versus the 
amplitude of the local random field h (/3 = 1, J = 1). The results were computed 
in the thermodynamic formalism using the new natural partition introduced by 
Behn and Lange ]l3| . The significance of the critical values hS"\ n = 1, . . . , 4, 
and /i£ 2a ' is explained in the text. The solid and dashed lines are exact lower and 



upper bounds on D q respectively which are obtained in section Q except for the 



d line coinciding with the values of D\ for h > which was obtained from 



/±( x ±) < \ an d converges to zero for f±(x±) > | leading to the critical value |ll| 

/j(3) = iarsinh(2-i(l -9e- 4/3J )5). (10) 

P 

The transition is again visible in numerically generated p n fl3| as well as in the D q - 
spectrum as the value of h for which D-oo begins to grow again for decreasing h, cf. 
figure 0. As was shown in [p9| , the generalized fractal dimension -D-oc has the value 

(3) 

Z?_oo = 1 at this point and _D_oo > 1 for h < he ■ 

The last of the transitions which are already well understood occurs at < 
when the slope of the coarse grained invariant measure density at x± jumps from q=oo 
to 0. The condition for this is f' a {x* a ) = Q resulting in |l3| 

h,W - iarsinh(3 • 2-i - \ - (3 • 2"i + \)e~^ J ) K (11) 

The transition is visible in numerically generated densities p n |l3| but not in the D q - 
spectrum, cf. hgure |[ Again, _D_oo can be calculated analytically and takes the value 
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Figure 3. Orbits of period three for different values of h. (a) shows the case 
where no point of the orbit falls into the overlap region whereas in (b) the point 
x*, ,,1 is in the overlap region and has therefore two predecessors (/3 = 1, J = 1). 



_D_oo = 2. In fact, Z?_oo can be calculated analytically for all h < he by considering 
that the scaling at the boundary in this case is weaker than at any other point ]To| . 
This gives lower and upper bounds on D q which for q — > — oo converge against each 
other, see section ||. 

Please note that the effects summarized in this subsection only depend on the 
measure at the boundary of its support and therefore are not strictly multifractal 
effects. 



2. Orbits and their contribution to the invariant measure 

The orbit to a given symbolic sequence {a} consists of all preimages frX ( x ^ a j), 
n G No- In the case of a periodic sequence {a} = ({<y}n)oa with finite period n, the 
orbit consists of the fixed points of the n functions f w i (T \ n in which tt denotes a cyclic 
permutation. To denote periodic orbits we will write for simplicity {er}„ instead of 
({ f7 }n)oo The fixed points x*i a \ of the periodic orbits are dense in the support of /i^ 
. Furthermore, any point of the support of /Zoo is contained in at least countably 
infinitely many orbits. 

2.1. Singularity of periodic orbits 

For the case of the fixed points x*± of f± it has been shown before that their singularity 
can be calculated explicitly [O, [l5| . Fixed points are 1 - orbits and we generalize this 
concept to periodic orbits of arbitrary period length. 

Let yi :— o o ■ • • o fa^( x \cr} ) , * = 1, . . - , n, be the points of the periodic 
n-orbit defined by {a} n . We then have y n — yo = x \a} because i^, is the fixed 
point of f{<r} n by definition. An example for a periodic 3-orbit is shown in figure |^. 

In case that no overlap exists or - if it exists - that no point j/i is in the overlap, the 
predecessor of each yi with respect to the iteration of the Frobenius-Perron equation 
is uniquely determined to be yi—\. There is therefore only one term in the Frobenius- 
Perron equation at all yi. We now investigate the singularity (pointwise dimension) 



W " + V W " " - * (12) 
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of /ioo at xj, . We assume| that the scaling limit 

- P °o( X {cr}„ + f ) _ ^(^{ct}* _ D 

exists for some finite k ^ and some ot{ a } n € As is shown in Appendix A.l| this 
implies 

, im ^(^J^K + f))-^>M.-t))^. (13) 

(CwJ'H.).) 1 ) " k 

The n-fold iteration of the Frobenius-Perron equation yields 

Poa(x*^ n + §) — Poo(a:| CT } n — f ) 

= H CT}n + D) - P-if^MrU ~ §)))■ ( 14 ) 

We denote the expression on the left hand side by X and the one on the right hand 
side by Y/2 n and thus have X/Y = 1/2". Inserting 1 = k/k and using (|l|) and (||) 
with the introduced notation X and Y we get 



1 



H [y X/e^^ J = 2"' ^ 

Most terms immediately cancel and we get 

((/wJ'(^wJ) QM " +1 = 2". (16) 

Therefore 

Tt In 2 

and with (/(-'(JVw.) = (/i„ } „ (/w.^.).)))" 1 = (/{.}. ^UX 1 " e «<* 

n In 2 

° w " = _1 " <*;.,.»• 1 ' 

This equation is invariant under cyclic permutation of {c}„ such that the scaling 
behaviour of the invariant measure at all yi, i = l,...,n, is given by the same 
Holder exponent a^ a j n . In other words, any point of the orbit has the same pointwise 
dimension D p = ct^ n —1. We therefore call a{ CT }„ the singularity of the orbit {cr} n . If 
a W\-n. < ^ *he measure has a positive (i.e. strong) singularity at all yi and if a{o-}„ > 
the singularity is negative (i.e. weak). 

As /wJ^jJ = n"=i/^(2/i-i)> the derivatives f^(y t -i) = A'(y 4 _i) at the 
points of the periodic orbit determine the singularity a{<x}„ and we finally have 

n In 2 

« M „=-l- ELilnA , (yi) - d9) 

§ Xt4 s assumption of strong scaling can be lifted in the non-overlapping case in which a generalization 
of (ha) to arbi trar y orbits can be proven provided the pointwise dimension of the fixed point exists, 
cf. subsection 2.2. 
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head (finite) tail (infinite) 

generic 

H h + H h- H h + H h- + H 1- ... 

W}n M 

periodic (here period 3) 

-+H H h+ - + H H M h + -++... 

{g} 9 = (- + +) 3 (- + +)oo 

offshoot (arbitrary head, periodic tail) 

H M h + - +H M M h + - + + ... 

{o-}„ (- + +)oo 

head determines the tail determines scaling at x*^ 
interval I{ a } n 3 x {a} 



Table 1. Illustration of the terminology of symbolic sequences and the different 
roles played by head and tail. 



A short calculation using the Frobenius-Perron equation (|5|) shows that the singulari- 
ties of arbitrary points x and fa-(x) are the same provided f a (x) is not in the overlap, 
cf. Appendix A.2| . The argument can be iterated such that the measure has the same 
singularity at any x and all its images f{a} m (x) for any {a} for which no point f{a-y i (x), 
i = 1, ... to is in O. Therefore not only the singularities at all points of a periodic orbit 
are the same but also the singularities at all points of non-periodic orbits in case the 
orbit also does not touch the overlap O. For orbits of the form {(T} m ({cr} I j) 00 we know 
that this singularity is the singularity of {cr}„, the periodic tail. We call non-periodic 
orbits of this type offshoots of the corresponding periodic orbit. The roles played by 
the head and the tail of a symbolic sequence {a} are summarized in table [l]. Note 
that the choice of the length of the head is arbitrary, in a sense. Similar structures 
have been considered in [E3|. 



2.2. Singularity of generic orbits in the non- overlapping case 

In the last subsection we have seen that the fact that if the predecessor of each point of 
an orbit is unique we can explicitly calculate the singularity of the orbit as a function 
of the derivative of A at the points of the orbit. In the non-overlapping case, O = 0, 
the uniqueness of predecessors holds for any point in supp /loo . It is therefore natural 
to try to extend ( |l9[ ) to generic non-periodic orbits. 
As a first step it is not difficult to show that if 

1 " 

lim -VlnA'(^ n) ) (20) 



n— >oc ri 

i=l 



exists for one choice of x ™ € f° r a given {a}, then it exists for any such choice 



and is independent of the particular choice made, cf. Appendix B 
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Let x G supp /Iqg be such that ( J20[ ) exists for the corresponding symbolic sequence 
{a}. To calculate the limit e — > in (||) it is sufficient to consider e n := |/{ CT }„| 
as 2|±i > min{A'(a;)|x G 1} = A 1 (x* + ) > holds. We know from the general 
properties that x — ij. G f° r au " G N. Because of the choice of e„ we 

have B En (x) D such that /j, 00 (B en ) > fJ.oc{I{a} n ) = leading to 

ln/ioo(gsj < -nln2 
lne„ ~ lne„ 

On the other hand we also can choose e' n := \f^ n (A)\. The interval I{ a } n is 
neighboured by two gaps. One of the neighbouring gaps is always /{ (T } n _ 1 (A), the other 
is either f{ a } m (A) with m < n — 1 or it is the complement of /. By contractivity of the 
RIFS we have \f{ a } m (A) | > |/{o-}„ (A) | for m < n, such that in cither case the smallest 
gap neighbouring 7 {o . }n is / {cr } n _ 1 (A). This implies ^(B^J < Moo(-f{a} n ) = 
because B e i can not bridge any of the neighbouring gaps and thus only intersects 
I{a} n - Therefore, 

In^ooCggJ -nln2 

lne; " h< ■ 1 j 

Using the mean value theorem for fs a \ n we obtain 

£n = |/ W J = /{,}„«) - / W >_) (23) 

n-l 

= (/W»)'(4 B> ) « - = II A ' (/{-}*«')) I'l (24) 
for some Xq™' G I. In the same fashion we get 

<=n^(/M 4 (*o cn) )) |a| (25) 

for some other a;^" 5 G I. Taking @ and (||) we obtain 

ln^oo^ej -ln2 
hmsup < lim sup- = (26) 

-nln2 

while using ( p2| ) and ( p5| ) yields 

In tioo (£L' ) — In 2 , „ 

lim inf , E " ; > lim inf , — 28 

™ " — iEiroln^(/ Wl K ( "»)) + iln|A| 

— «ln 2 . . 

lim — -. = r . (29) 



n — >oc 



In both ([27]) and (^9|) the existence of the limit in (|2^) and the independence of ( pp| ) of 
the points a;^ G was used to replace lim inf and lim sup by lim and to substitute 
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the points of the orbit {a} for Xq" and x' (rl) respectively. From (^) and ( p9j ) we 
immediately get that D p exists and is given by 

— nln2 

D n = hm — = ; — r- (30) 

Elton's ergodic theorem (3^] implies that our assumption of the existence of (pcj) holds 
for almost all {cr} which corresponds to fioc - almost sure existence of the pointwisc 
dimension D p (x). Elton's theorem together with ( |30| ) further implies that for pb^ - 
almost all x the pointwise dimension takes the common value 

^ioo-a.s. — In 2 

D pW = rTl ^777^ 7ms= :D p- ( 31 ) 

J (In^'(O)Moo(^) 

This result has direct consequences for the information dimension D\. Proposition 
2.1. in |3(J] implies that 

J} ff ({ieJ:flp(i) = flp})='D p (32) 

in which Djj denotes the Hausdorff-dimcnsion. General properties of the multifractal 
/(a)-spectrum imply that a% is the only fixed point of f(a) and that a% — D±, cf. 
$7\ . Therefore ([52j) implies £>i = Z3 p . This fact is illustrated by the solid line in 
figure U coinciding with the numerically obtained values of D\ for h > he which 
was obtained through calculation of ( |3l| ) using Edalat's i?-integration S). (It also is 
a simple exercise to check on a computer that using the first 10 5 digits of the dual 
representation of it or e as a symbolic sequence, the value of a obtained from ( |30| ) 
is also exactly the value D\ — 1 in figure ||. Of course, the use of a random number 
generator instead of 7r or e yields the same result with probability one.) 

Please note that the restriction to almost all {cr} in the above is necessary as the 
sum in ( ]20| ) does not converge for all {a}. A simple example in which it does not 
converge is a sequence of bulks of plus and minus signs of ever increasing length. The 
length of the bulks can be chosen such that the sum in (12(3) keeps oscillating for any 
size of n. 



2. 3. Singularity of orbits in the overlapping case 

In the previous subsections the condition that no point of the orbits under considera- 
tion is in the overlap O was essential for the calculation of their singularity. In this 
section we investigate how the singularity of orbits is affected by the overlap O. Tuning 
the parameter h changes O as well as the location of the orbits. If (at least) one point x 
of a periodic orbit is in O this point has two predecessors and there are two terms in the 
Frobenius-Perron equation (j^) at this x. The two terms contribute singularities a± and 
a 2 , i.e. m := fi 00 (B £ (fZ 1 {x))) ~ e ai and fi 2 ■= Moo(S e (/+ 1 (a;))) - e° 2 . Therefore, 
the singularity at x will be inf{ai,a2}, since fj, 0o (B e (x)) ^^1+^2 ~ £ aJ t a ^ a '}. 
The mechanism is illustrated in figure |i| using the example of an orbit which will be 
important in the next subsection. In the case that the original singularity is stronger 
than or equal to the additionally contributed one, there are no consequences. In the 
case that the singularity is rather weak though, a weak singularity is replaced by a 
stronger singularity. In fact, the new singularity at the maximal value of h for which 
x is in O is always rather strong as it stems from x± where A' is small. The change 
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iteration depth 



Figure 4. This figure illustrates the intricate structure of orbits and how weak 
singularities are superseded by stronger ones. The point x = x* s , , , is 

f2) (2a) y°°i 

in the overlap at /? = 1, J = 1, h = 0.74 (h c < h < h c ). In this figure x, 
all its predecessors up to 6-fold application of and all its successors up to 
2-fold application of f± are shown. The points are the predecessors/successors 
whereas the connecting lines illustrate which point is mapped onto which. Note 
that the backward trajectory of x branches each time a point is in the overlap 
whereas the predecessor is unique if the overlap is not touched. The dashed line 

connects points of the orbit {H (H — )ac} which is an offshoot of the {-I — } orbit 

and therefore has the singularity j « 0.982. The solid lines connect points 

of two other orbits which meet at y G O, both carrying the generic singularity 

a = Dp — 1 f» —0.057. At x the weaker singularity j is superseded by a. 

The offshoots emerging from x all have the stronger singularity a. 



in the singularity may have a major impact on the -D^-spectrum especially if some or 
all of the weak but somewhat stronger singularities have already vanished. 

A special role in the mechanism described above is played by the 1-orbits {+} 
and {— }, because they never touch O, and the 2-orbit {-I — }. At moderate overlap 

(small \0\) we have the situation illustrated in figure ||. Since x*^ j is mapped to 

x*^ |_i and vice versa and because f a is monotone, all points to the right of x*^ j are 

mapped to the right of x| ^ and all points to left of x*r ^ are mapped to points 



left of x*^ + _ 



y Therefore, any periodic n-orbit with n > 2 must have at least one point 
J. Hence, the {-I — } orbit is the last periodic orbit to be reached 



"{-+}> x {+ 



byO. 

The {-I — } orbit and its offshoots carry a very weak singularity as the {-I — } orbit 
always stays in regions with comparably large A 1 , cf. equation (19). For (3 and J in 
the vicinity of f3 = J = 1 one can show that the orbit {H — } has even the weakest 
singularity of all periodic orbits . Because of its weak singularity and the fact that 

all other periodic orbits and their offshoots are reached by O before x*^ j is reached, 

the {-I — } orbit and its offshoots practically solely determine all D q with q < if h 
is such that O has nearly reached x*^ y If O includes x*^ y i.e. /_(a;^_) > x*^ 

" J ■ ' i) the preimage /r 1 (a;| + _j). (The 
of course.) The additional preimages contribute a stronger 



x'^ -j has additionally to the preimage / + (x^ + 



same applies to x 



{-+}' 
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/- /+ 

Figure 5. Mapping of subintervals of / = [x*_, x*A under /_)_ and /_ elucidating 

the importance of the interval [x^ ,]. The latter is mapped onto parts 

of the outer intervals [a;?l,a;| ^j] and ji^+l under /_ and /+ while these 

outer intervals are itself preimages of parts of [x^ |_, }!■ ^ ne P om ts in 

the overlap f—(x^_)] have two predecessors, one stemming from the left 

under /_)_ and one from the right under /_ (all in the case of moderate overlap). 



singularity than the original weak singularity of the {H — } orbit. Thus, the weak 
singularity of the {H — } orbit and its offshoots is superseded and, as all other periodic 
orbits have been reached by O before, all D q with negative q have collapsed to D q = 1 

(2) 

at this point. The critical value h c at which the collapse takes place is therefore 
given by the condition 

f-(xl) = x* {+ _ } . (33) 

So far we have only discussed periodic orbits and their offshoots. Other non- 
periodic orbits do not play a major role because they generically have the rather strong 

singularity a — D p — 1 = D\ — 1 and also generically have points inside [a;| ^ , x*,^ J 

such that they are reached by O before the {H — } orbit is reached. 

After the collapse the right part of the /(a)-spectrum of the invariant measure 
has vanished since the weaker negative singularities have all been superseded by 
stronger ones, cf. figure |6^. A similar collapse of parts of the multifractal spectrum 
has previously been observed in the superposition of equal-scale |59| and multi-scale 
pp| , |4l|] Cantor sets showing that effects of this type appear in a wide variety of 
applications. 

(2) 

Let us now determine h c explicitly. In a first step we need explicit expressions 
for x+ — —x*_ and x|_j , = — x*t ^. The fixed point x+ is defined by f+(x*^_) = x* + . 

With the notation z = e 2 @ x + this yields the equation 

z 2 - (e^ J (e 2l3h - l))z - e 2f3h = (34) 



with the solution 



= -!- In (K + VK 2 + e 2 ? h ) (35) 
2p 




Figure 6. Collapse of the right part of the /(ct)-spectrum at and h^ a \ 
In (a) the spectra at h = 0.49385 > h {2) (dashed line) and h = 0.4938 < 
(solid line) are shown. In (b) the spectra are at h = 0.8136 > h^ a ^ (dashed line) 
and h = 0.8128 < h{ 2a) (solid line). The spectra were obtained by a numerical 
Legendre transform of the corresponding D^-spectra. These were generated with 
the same algorithm as in figure fcj with a recursion depth of 21 (/3 = 1, J = 1). 



where K — e 2/3J ' (e 2 P h — l)/2. To obtain we exploit 

-f-(x} , _■>). With z = e 2Px i+-> this yields 



~ x {-+) 



and therefore 



z 2 + e -20J {e 20h_ 1)z _ e Wh = Q (36) 



x* {+ _ } = ^ln (K + VK 2 + e^ h ) (37) 

with K = e -2/3J( e 2/3h _i)/2. 

The equation f-(x* + ) = x*,^ j is equivalent to x* + — x*,^ j + 2h. With the 

explicit expressions for x* + and x*^ we get 

2(cosh(2/3^)) 3 + 3(cosh(2/3/i^)) 2 = ( cosh(2/3J)) 2 . (38) 

This equation has exactly one real solution for cosh (2f3h^) resulting in 

/i< 2) = — arcosh( cosh(§/3 J) - §) . (39) 
zp 

The phase diagram for the transition at is shown in figure 0. One clearly sees 
that there is a critical line /3( J) such that there is no phase transition possible for all 

(2) 

(3 < j3{J). This line is given by the condition h c = corresponding to 

P(J) = ^arcosh(|) = -Lin (2 + V5). (40) 
With (^9|) and ( |40| ) we thus have a complete analytical understanding of the occurrence 

(2) 

and the position of the transition at he ■ 
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,(2) 



(2) 

Figure 7. Surface plot of the critical value h). ' against /3 and J. The plot was 
obtained from 



2.4- Transition of the measure at h^ a ^ 

We have seen in the last subsection that the sharp drop of all D q with negative q at 

h c can be explained through the analysis of the properties of the {H — } orbit which 
is on the one hand an orbit with very weak singularity and on the other hand the 
last periodic orbit (with decreasing h) to be affected by th overlap O. To understand 
the smaller drop at h c 2a ^ it is therefore natural to look for periodic orbits with weak 
singularity which also are affected by the overlap relatively late. We will now argue 

that the family of periodic orbits of the form {H (H — )„}, n £ N fulfills these 

conditions and that their change in singularity indeed causes the transition at h^ a \ 
To present the precise argument we need to investigate some properties of this family 
of orbits. 

For large n most points of the {H (+— )„} orbits are close to the points of the 

{H — } orbit. Thus, having ( |l9| ) in mind, the singularities of these orbits have a similar 
strength as singularity of the {H — } orbit, i.e. they are weak, provided that no point 
of the orbits is in the overlap O. The larger n the larger is the fraction of points of the 

corresponding orbit which are close to x*^ j and x*^ ^. Therefore, the singularities 

get weaker with growing n. 

Concerning the position of the orbits {H 

z n := x*j ^ j j are for each of the orbits {- 

As / {+ — (+-)„}«') 



J {+- 
overlap O 



„(») 



T+— (+-)< 



^ )«} we first note that the points 
-(+—)„} the closest points to the 

j for any choice of initial points 



G / and n — > oo, we get with Xq := z n = 



that 



'{+-(+-)») 



for n — * oo. Furthermore, it is an easy exercise to check that z n is monotonously 
growing with n. 

Taking these properties together we conclude that each of the periodic orbits of 
the form {H (H — )„} is affected by the overlap as soon as z n is in O and that the 
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1.1 1.2 1.3 f3 



Figure 8. h c obtained from the analytical condition ( pjj ) (dashed line) 
compared with the location of the small drops in D q (error-bars), cf. figure 0. 
The errors are estimated and mainly due to the errors in the determination of the 
location of the drops in the numerically obtained D 9 -spectra (J = 1). 



orbits with weaker singularity (those with large n) are affected by the overlap later 
than the orbits with somewhat stronger singularity. Finally, as the z n converge to 

x*^, ^ ^ j, there are still countably infinitely many orbits {H (H — )„}, n > N 

with some N € N left which are not affected by the overlap as long as x^ ^ y is 

not in O. As soon as x^ + __^ j , gets into the overlap O, all orbits {H (H — )„} 

are in the overlap and their singularity therefore superseded. The critical field strength 
is thus determined by the condition 

/-«) = <-(+-)„ = f{+-}( x {-+})- ( 41 ) 

This criterion is in perfect agreement with the numerically obtained positions of the 
small drop in the D g -spectrum for q = —1, —2, —3, cf. figure ||. The transition is also 
visible in the /(a)-spectrum. For h — > h^c°^ + a cusp develops and the spectrum 
collapses to a smooth form again at h^? a \ cf. figure ^b. 

So far we only addressed the periodic orbits {H (H — )„} neglecting their 

offshoots. Each such orbit has countably many offshoots. From these offshoots those 
originating from points left of O and containing exclusively additional — as well as 
those originating from points to the right of O and only containing additional + are 
also not affected by the overlap as long as the corresponding periodic orbit is not. 
Thus, at ft.i 2a ' countably infinitely many periodic orbits with each countably many 
offshoots of the described form vanish at once. This explains why this transition is 
visible in the Z) g -spectrum in contrast to events of single orbits (and their offshoots) 
being affected by the overlap at any value of h. 

The fact that the transition is not visible in D q with large negative q is also easily 
understood in this framework. As the singularities of all {H (H — )„} orbits are 
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Figure 9. Double logarithmic plot of p n generated from the band structure after 
n = 21 iterations. The negative peaks correspond to deep cuts in p n . The dotted 
lines mark the position of iE|f_|_w_|_ \ i, i.e. the points of the {H — } orbit and 

its offshoots, and the dash-dotted lines the positions of x* rl , . . , . , , , i.e. 

U+JaH (+—)oci' 

the points of the {H (-) — )oo} orbit and its offshoots. In (a) for h > h^? a ^ 

deep cuts are visible at all marked positions whereas in (b) for h < /ic 2 "' the deep 
cuts at (+— ) } nave vanished. They reappear in (c) at a critical field 

strength h ss 0.5436 as is discussed in the text whereas in (d) for h < he all deep 
cuts have vanished Q3 = 1, J = 1). 



similar to but somewhat stronger than the singularity of the {H — } orbit, the {H — } 
orbit and its offshoots dominate D q for large negative q and the transition is not 
visible. 

Please note the symmetry of the system which allows the same reasoning with 
the 'opposite' orbits {— + + ( — h) n } resulting in an equivalent result. 

Even though the effect is due to the periodic orbits {H (+— )„} and their 

offshoots, the orbit entering into condition (|4l| ) is an offshoot of the {H — } orbit, 
namely {H (+— )oo}- This orbit thus seems to play a similar role as the {H — } 

(2) 

orbit for the transition at he . For illustration we have generated an approximation 
of p n and calculated the position of the points of the { +— } orbit and its offshoots of 
the form {(+) n (H — )oo} (dotted lines in figure ^|) as well as the position of the points 

of the {H (+— )oo} orbit and its offshoots of the form {(+) n H (H — )oo} (dashed 

lines in figure ^|) . It is obvious that there is a bunch of weak singularities in the vicinity 
of x|j f, ■) l and in the vicinity of the points of the offshoots which all vanish when 
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Figure 10. Phase diagram for the transitions of the invariant measure derived 

(n) 

from the random field Ising model. The h c ' are the critical random field strengths 

(2) (3) 

defined in the text. The remarkable fact that the lines of h c ' and h\ 1 as well 
as the lines of he and /i' 3 ' intersect shows that there is a variety of different 
scenarios depending on the choice of kT. For example at kt = f .3, the transition 
at he precedes the one at h c 3 ^ while the transition at h c 4 ^ is non-existent ( J = f ). 



decreasing h below h c 2a \ cf. figure [)]a and b. 

The attentive reader will have noticed that we have argued that the weak 
singularity of an orbit which is touched by the overlap O for some ho but not for h > ho 
will be superseded by the strong singularity of the {+} or the {— } orbit. For smaller 
h in the generic case the additionally contributed singularity will be D\ — 1 and thus 
also rather large. It is not excluded though that two weak singularities are combined 
resulting in a weak singularity even though the orbit is in the overlap. A prominent 

example of this effect is the value of h for which x|_| ^ \ i = x*< ^ , = 0. 

At this value of h new deep cuts in the approximated measure density appear at 

x*^ (+-)«,} ana - offshoots, cf. figure §c. The effect on the Dq-spectrum is 

negligible though as this is a rare event only affecting a single orbit and its offshoots 
at a given h. 

To summarize, the crucial feature for a visible transition in the Z) g -spectrum is 
that a non-negligible fraction of orbits with weak singularities is affected by the overlap 
at one sharp critical value h c resulting in a drop of D q with negative q. It is not 
excluded that there are more transitions of this type which might become observable 
with further increase in numerical accuracy in the future. Our arguments hint to the 
conjecture that these transitions, if existent, should take place at ft, > h c 2a \ 

The formulae given in the last two subsections allow to draw a phase diagram for 
the phase transitions in the -Dg-spectrum of the invariant measure which have been 
observed so far, cf. figure [l^. It should be emphasized that these 'phase transitions' 
are not phase transitions of physical quantities like the (local) magnetisation or the 
Edwards- Anderson parameter. 
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3. Bounds on the .Dg-spectrum 

Let in the following /i denote the invariant measure. As mentioned above there are 
natural bounds on D q . They are induced by 

/^)<£/4 (42) 

i 

where x is some point in the support of fi and n^ x ) is the measure of the box containing 
x. Let us now first consider the case q < 0. In this case inequality ([l2]) immediately 
yields 



£> g > r ^- r lim ™ J . (43) 



q lim MfJjz 

q — 1 e^o In £ 

The limit on the right hand side is the pointwise dimension at x. If we choose x = x* + 
we can calculate the pointwise dimension at x for any h,T ^ 1, J using (|l^). This 
gives the lower bound 

q In 2 
l^gln(,4'(a;* ))' 



D i > i ( 44 ) 



D i ^ T^ZTZTTTTZT* W ( 4 5) 



(2) 

In the region h > h c we also can calculate the pointwise dimension at x*^ i. With 

x = | we get the bound 

In 2 

These lower bounds are shown as solid lines in figure 0. 

Whenever we know that the invariant measure scales most weakly at some point 
x of its support and we can calculate the pointwise dimension at x, we also can give 
an upper bound induced by 

E^M-v ( 46 ) 

i 

In this N is the number of boxes and is essentially proportional to 1/e. In the region 
h < hf ] the measure is assumed to scale most weakly at x* + (and x*_) such that 
inserting ( [i"6|) into (Q) we obtain the upper bound 



1- ? V ln(A'(i-;» 

(2) 

cgion h > he the scalir 

yielding 



(2) 

In the region h > h c the scaling is assumed to be weakest at a? |_| j (and air |i) 



+ <48 » 

Obviously the upper and lower bounds converge against each other in both regions as 
g — > —00 such that we get the explicit expressions 
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The expression for the region h < was already given in p9| whereas the expression 

(2) 

for the region h > h y c ' needs our analysis of the singularity of orbits in section 



2.1 



These upper bounds are in general not as good as the lower ones found above and are 
not shown in figure ^. 

In the case q > all considerations remain valid for the bounds based on ( ff2| ) 
which yields by use of ([l9|) the upper bounds 

q In 2 

D q < 50 
1 -q]n(A'(x%)) 

for any h and the sharper condition 

q In 2 

D " ~ wma'^j)) (51) 

for h > h [ c ] . These upper bounds © are shown as dashed lines in figure 2. 

Lower bounds on D q can be obtained by considering points x at which the scaling 

(3) 

is maximally strong. This is assumed to be true for x* + and x*_ in the region h > h c . 
Therefore we get 



gln2 

1 - q V ' ln(A(x~%)) 



P«>— (1 ( 52 ) 



These lower bounds are not as good as the upper bounds ( pl| ) and are not shown in 
figure ||. In the region h > the lower and upper bounds for D q converge against 
a common value for q — > oo yielding 

In 2 

D °° = ~ha>(x* + )Y (53) 

The expressions (^) and ([33|) generalize results priorly obtained [|l0| for h > . 
Note that for q < (g > 0) the upper (lower) bounds rest on the assumption of 

minimal (maximal) scaling at x* + {x*_) and x*,^ j (x| |_i) whereas the lower (upper) 

bounds do not need any additional assumptions. In figure ^ only the lower (upper) 
bounds are shown. 

The fact that for q < — 2 (q > 4) the lower (upper) bounds are more or less 
identical with the numerical data (cf . figure ||) shows that in these cases the generalized 
dimensions D q are solely determined by the scaling at certain points, i.e. at x* + and 

x*_ for h < h^ and at x*^ , for h > h^K This reinforces our previous argument 

why the small drop in the £> g -spectrum at h 1 ^^ is not visible in numerical data for 
q<-3. 



4. Concluding remarks 

By generalizing arguments, hitherto mainly applied to fixed points, to orbits we have 
been able to calculate the singularity of all periodic orbits not touching the overlap 
O and the generic (/i-a.s.) singularity of orbits in the non-overlapping case. We then 
investigated the effects of a non-void overlap on the singularity of orbits. While being 
a relevant result in its own right the knowledge of the singularity of orbits and their 
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Figure 11. Sample situation of first order bands in the symmetric four Dirac 
mass case in the vicinity of a transition similar to the one at he ' in the text. If, 
as shown in the figure, H is large enough, £| 2 3} an< ^ x {32} are no * mc l u ded in 
I\ and I2 respectively. In this case £{ 2 3} nas * ne umc l ue predecessor a;| 32 ^ and 
vice versa exactly like x*^_ + y and in the dichotomous case. The situation 

is therefore identical to the situation in the dichotomous case and the transition 
takes place at h = ■ (/3 = 1, h = 0.7, H = 1.8) 



dependence on touching the overlap or not also provided the explanation for two phase 
transitions in the multifractal D^-spectrum of the invariant measure. 

RIFS with similar properties as the one discussed here also appear in a variety of 
other contexts [pTj, ^2], |2^, 29 1. In |29j a similar phase transition in the Z? 9 -spectrum 



already has been observed. The explanation of the mechanism causing such transitions 
does not crucially depend on the special choice of the RIFS we investigated here. Let 
us briefly discuss some obvious generalizations. 

Dropping the restriction of symmetry of the random field distribution leads to the 

loss of various symmetries of the orbit structure, e.g. x* + 7^ — x*_ and x*^ , ^ —act h 1 , 

and of the resulting calculational simplicity. The general orbit structure underlying 
our explanations of transitions in the D 9 -spectrum however remains. Therefore, the 
transitions also exist in the non-symmetric case and can be analysed in the same way 
as above in the symmetric case. Some numerical results may be found |jl3| . The 
explicit calculations are far more complicated though. 

Considering more complex random field distributions like e.g. n (symmetrically 
or non-symmetrically distributed) Dirac masses considerably complicates the system. 
The mechanisms causing phase transitions in the D 9 -spectra are nevertheless still 
the same. Let us briefly outline the example of four symmetrically distributed Dirac 
masses located at {—H,—h,h,H} as the random field distribution, i.e. the RIFS 

h, h, fi} = {A- H,A-h,A + h,A + H} with h < H € K+: 

• For any given [3 a transition line in the /i-7?-plane exists between a region with 
Dq = 1 and Dq < 1 which is easily determined from the overlap conditions 
of the first order bands Ij, j = 1, . . . ,4. These conditions read h < and 
h{ x X) — fi( x i)t t ne second inequality being a condition on H as well as h. 

• The transitions depending on the scaling at the boundary of the support of /Xoo , 

i. e. on the scaling at x\ and exist in the same way as those at and 
above. The transition conditions are H = h^P and H = respectively. 

• Even though the orbit structure is more intricate than in the dichotomous case 
figure pi] shows a situation in which a transition of the type of the one at he 

— (2) 

above can take place. The transition occurs at h = hi provided H is large 
enough for the bands Ji and 1% not to include a;| 2 3} or ^{32}' ^- e - fi( x t) < x {23} 
or equivalently fi{x\) > x*^ 32 y 
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In summary, the types of transitions occurring in the D g -spectrum are the same as in 
the dichotomous case and it is possible to draw a diagram in the {(3, h, H) parameter 
space similar to the simpler diagram in figure [l^ which is the h = H slice of this higher 
dimensional diagram. 

In the same fashion the analytic tools developed in this paper in principal allow 
to analyse the Z? 9 -spectrum of the distribution of the effective field of the RFIM for 
any discrete random field distribution. Moreover, the exact form of the function A is 
not crucial such that the analysis can be performed as soon as the following features 
are present: 

• A hyperbolic RIFS {/?}, j — 1, ... ,n which - in dependence on a control 
parameter h - has overlapping bands or not. 

• The functions fj are sufficiently smooth, monotonous and have comparably large 
derivatives at certain periodic orbits. 

• The conditions for the transitions that certain (periodic) orbits touch or do not 
touch the overlap can be fulfilled by tuning h. 

The generalization to continuous random field distributions is not obvious as the 
techniques used above can not directly be applied to this case. In IJ a numerical survey 
of the nature of the distribution function P^ for various random field distributions 
can be found. For discrete random fields a transition between a devils's staircase 
corresponding to Dq < 1 and a smooth function corresponding to D$ — 1 is found as 
is predicted analytically. For continuous random field distributions without gaps P^ 
always is smooth whereas gaps in the continuous random field distribution define a 
scale above which (for certain parameters) P^ resembles a devil's staircase and below 
which it is always smooth. 

The generalization to more general lattices like the dichotomous symmetric RFIM 
on the Bethe lattice (cf. p2| E3f and references therein) is non-trivial. Transitions 
corresponding to those at h^^ and he above are present. The exact conditions 
determining the critical field strengths are however not obvious and need further 
careful investigation. The existence and properties of transitions of the type of those 
at or hi? a ^ needs even more careful and detailed analysis and must be deferred 
to further work. We stress however that all these transitions in the Z^-spectrum take 
place far in the contracting, paramagnetic regime of the physical phase diagram of the 
RFIM on the Bethe lattice. 

The bounds on the D g -spectrum obtained in section |^ are also of a quite general 
type and therefore applicable in a variety of contexts. As figure || shows, the 
bounds together with the explanation of all the transitions discussed give a very good 
qualitative as well as quantitative understanding of the D g -spectrum of the invariant 
measure of the effective field in the dichotomous, symmetric ID RFIM. 

The detailed analysis of the invariant measure of the effective field should be 
viewed as the preparation of the study of the multifractal properties of the measure of 
the local magnetisation in the RFIM which essentially is a convolution of the invariant 
measure of the effective field with a distorted version of itself. The transitions in the 
Z? g -spectrum of the invariant measure of the effective field have direct counterparts 
in the Dq-spectrum of the measure of the local magnetisation and are thus gaining a 
direct physical significance. We will address this subject in forthcoming work. 
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Appendix A. Relations between scaling properties 

Appendix A.l. Relation between the scaling at x*^ and at f7 a \ ( x {a} ) 
In this Appendix we show that the assumed scaling relation 

lim P ^ g W- + s)-y g W--§) =: k (AJ) 

implies that the limit 

lim ^(/^^)„+f))- f (/M„^k-f)) (A 2) 



exists and is equal to k. In the following we drop the indices {<r} n and oo to improve 
readability and denote the expression in (A. 2) as Q(e). Applying the mean value 
theorem to / _1 in Q(e) and using f~ 1 (x*) = x* yields 

OM p (-* + (TW + m p(x* (f-iy(x« - fr)f ) 

with some 5i,5 2 E [0, |]. Now let (J -1 )^ De tne minimum of (}~ r )'{x* + 6i) and 
— 5%). Because / _1 is strictly monotonously growing we have (/ _1 )min > 0- 
Using the fact that P is monotonously growing as well we get the lower estimate 

n(A > ^ + (r 1 Linf)-^-(r 1 Linf) rA 

((/ 1 )'(!!!*) e) 



The quotient of (/ -1 )mm an d converges to 1 and thus 



(A.5) 



lim Q{e) > lim ^^^^l^fi = *• (A.6) 

Using the maximum of (/ _1 )'(a;* + Si) and (/ _1 )'(a;* — #2) instead of the minimum 
we get in the same fashion the upper estimate 

lim Q(e) < k. (A.7) 

6— >00 

Both estimates together give the conjectured result lim e ^ Q(z) = k. 
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Appendix A. 2. Relation between a(x) and a(f~ (x)) 

The Frobenius-Perron equation for the invariant distribution induces the equality 
a(x) — ct(f~ (x)) for any x G supp/i which is not in the overlap O. 

The proof is straightforward. Let x £ supp /i be a point which is not in O and 
a G { — , +} an arbitrary sign. Then, dropping again the index oo, 

a(x) = hm - — (A. 8) 

e^a lne 



1S _ H±P(f-\x + §)) - \P{f-\x- §))) 



e^O In e- 



(A.9) 



because of the Frobenius-Perron equation. We now use the mean value theorem for 
+ |) and f^ 1 (x — |) to obtain 

r -In 2 , HP(f-^x) + (f-^(x + S 1 )^-P(f-\x)^(f-^(x^6 2 )^)) 

= hm — 1 — 

e^o In e In e 

(A.10) 

with some <5i , ^2 € [0, §]. Using again the notation (/,7 1 )min f° r the the niinimum of 
{fa l )'i x + S i) and (/a 1 )'^ - § z) and defining e' := (Z^ 1 )^ • £ we get the inequality 



HP(f- 1 (x) + ^)-P(f-\x ) - 
s'-0 -ln(/ CT - 1 ); nin + ln £ ' 



Using the maximum (/ CT 1 )J nax of the derivatives instead of the minimum we get by 
the same token a(x) < a(f~ 1 (x)). Therefore equality follows. 



Appendix B. Independence of (20) from the choice of a^™' 
In this Appendix we show that 

1 n n 

hm ~(£]n0OO-5>#5i n) )) =0 (B.l) 

n— too n \ — '—^ ' 

i—1 i—1 

for any strictly positive diffcrentiable function </>, any symbolic sequence {a} and any 
choice of x-"',^™' £ Ifa}^ i 6 {1, • • • n £ N. Let 0, {c}, x-"' and x-" 1 be given 
and set e n :— \I{cr} n I- Then, using the mean value theorem, 

0(5?-') = <M*n + w° + - ^) (b.2) 

with some 5™, |(5-"'| < Thus 

ln0(a*») - ln0(x<"') = In (l + ^-±|7l ( ^) _ 8 W)) . (B .3) 

We denote the finite constant max { : x € /} by Q max . For sufficiently large 
i £ N the expression 1 — Q max |x^' — x-"'| is positive and we get 

hi (l - Q max |xf -a£°|) < ln0(x<' l) )-ln0(x<' l, ) < In (l + Q max |x<"> - x<"> |) . 

(B.4) 
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With \x l " } — a;-"'! < er, this yields 

ln(l - Q max £i) < In <p(x\ n) ) - ln^(a£°) < ln(l + Q max£l ) (B.5) 
implying that the difference \n<f>{x^) — hi(j)(x^) converges to zero. A standard 



argument then shows that the average also converges to zero, i.e. that (B.l) is true. 
With cf> — A' this yields the alleged independence of ( pp| ) from the choice of x^ . 
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